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A general scheme for rotational cooling of diatomic heteronuclear molecules is proposed. It uses a 
superconducting microwave cavity to enhance the spontaneous decay via Purcell effect. Rotational 
cooling can be induced by sequentially tuning each rotational transition to cavity resonance, starting 
£f) ■ from the highest transition level to the lowest using an electric field. Electrostatic multipoles can be 

used to provide large confinement volume with essentially homogeneous background electric field. 
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I. INTRODUCTION 



Currently, there are three main techniques to produce translationally cold molecules. The buffer gas [1] and Stark 
deceleration [2] cooling schemes employ trapping of the low field seekers which requires the molecules to be in the 
! , appropriate internal states. Cold molecules produced from photoassociation technique [4] are usually vibrationally hot, 
' occupying a number of high lying vibrational states. Vibrational cooling schemes using optimal control of ultrashort 
^sO i pulses have been proposed [5]. Translationally cold molecules are useful for precision spectroscopy and measurements, 
molecular optics and interferometry as well as cold collisions studies, although they may be internally hot. It is also 
desirable to have translationally and internally cold molecules as well if the aim is to obtain molecular Bose-Einstein 
Condensate [3]. For the purpose of trapping, the molecules need to be cooled internally. 

There has been no scheme to rotationally cool molecules although there are several methods of producing molecules 
in a single rotational state. A classic example is the electrostatic low-J selector which has been used for rotational 
state selection long ago [6] . Optical Stern-Gerlach has been demonstrated for atoms [7] and molecular state selection 
by lasers has been proposed [8]. However, these techniques are best applied to molecular beam. They are not cooling 
schemes since they do not employ a dissipative mechanism. Supersonic expansion is based on dissipation and may 
also produce rotationally cold molecules, but cannot be used to cool molecules which are already translationally cold 
since it requires high pressure and can only be applied mechanically on hot molecules. Besides, for translationally 
cold molecules it is more efficient to employ dissipative process for internal cooling which avoids the loss of molecules. 

Recently, we have proposed 1-D translational cooling schemes for molecules which rely on a single optical sponta- 
neous emission [9] . The schemes can be repeated for 3-D cooling if the 1-D cooled molecules can be brought back to 
the initial internal state again. Molecules in excited rot-vibrational states can redistribute significant amount of the 
internal energy into the translation motion by state-changing inelastic collisions. At high temperature, this becomes a 
problem for trapping due to the increased collision rate. Therefore, it is essential to removed the internal excitations 
by spontaneous emissions fast enough before the inelastic collisions occur. Spontaneous emission carries away entropy 
to radiation while inelastic collision transfers entropy to the translational degree. For sufficiently dilute heteronuclear 
molecules, vibrational cooling occur within 0.1s — Is timescale through infrared spontaneous emissions before the 
vibrational inelastic collisions take place. Thus, the vibrational entropy is discarded as radiation entropy instead 
of the translational entropy. However, rotational spontaneous emissions can occur only in polar molecules and take 
much longer time, beyond the experimental timescale for dilute gas. For dense gas, rotational decay may be enhanced 
through many-body effect, but inelastic collision rate may become dominantly large. 

In this paper, we propose a rotational cooling scheme for confined polar molecules. We use external electric field 
to tune the internal transitions into resonance with a lossy microwave cavity (Fig. 1). This enhances the rotational 
spontaneous emissions via the Purcell effect [10] whereby the populations are transferred to a lower rotational level 
in stages and eventually to the ground rotational level. The Purcell effect has been experimentally demonstrated by 
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Goy et. al. using Rydberg atoms in a microwave cavity [11]. The cooling process requires long interaction time and 
we have proposed a method to confine the molecules within the cavity field. 

II. INITIAL STATE 

We consider a gas of polar diatomic molecules at thermal equilibrium temperature T with the Maxwcll-Boltzmann 
distribution. We also assume that the gas is sufficiently dilute such that the inelastic collision is negligible throughout 
the cooling process. At below room temperature, the populations are essentially in the eigenenergy in the ground 
electronic state \g,Q). Thus, the population and the eigenenergy in the vibrational state \n) and rotational states 
| J, M) with J > Q (the total electronic angular momentum) are given respectively by 

P(g,n,n, J) = P (2J + 1) ]T evp{-E g ,n,nj/k B T} (1) 

2 ! 2" 

E g , n j,n = hc{T g + co g (n + 1) - w gXg (n + \f + (B g - a g (n + ±))(J(J + 1) - O 2 )} (2) 

where T g ,uj gi uj g x gi B g and a g are the standard electronic, vibrational and rotational constants (in wavenumber) [12]. 

From Eqs. 1 and 2, we estimate that for CsF molecule with (u> e , aj e x e , B e , a e ) = (352. 56, 1.61, 0.184, 0.0012)cm _1 
[12], only the first 5 rotational levels in the ground vibrational level [n = 0) arc significantly occupied at temperature 
around IK. The rotational levels and the magnetic (Zeeman) states for cases O = and fi = | are shown in Fig. 
2a and Fig. 3a respectively. Initially, the molecules are in all the magnetic states. In order to cool the molecules in 
all the magnetic states, the cavity must be able to support linearly polarized as well as circularly polarized photons 
corresponding to AM = and ±1 transitions respectively. If we assume that the cavity supports only linearly 
polarized photons, only the states with M < ±£l will be sequentially Stark-tuned to cavity frequency by an external 
static electric field for enhanced spontaneous decay toward the ground level. Molecules in other states (with M > ±0) 
are off-resonance with the cavity and their decays will not be enhanced. In this case, we will have to extract only the 
molecules in the states M < ±f2 for rotational cooling. 

III. ROTATIONAL COOLING SCHEME 

In this section, we elaborate on the rotational cooling scheme based on the proposed setup apparatus (Figure 1). 

A. Cooling mechanism 

The spacings of rotational levels in polar diatomic molecules increase monotonically with energy levels. We use this 
monotonic property to bring the populations from the high J levels to a low J (ground) level by enhanced rotational 
decay in a lossy microwave cavity. We want to bring each pair of transition (| J, M) to | J — 1, M — q) with q = 0, ±1 
and J ^ 1) into resonance with the cavity at each time, starting from the J max — » J max — 1 transition and ending 
with the J = fl + l^J = £l transition. We shall elaborate on the simplest scheme employing only the 7r-transitions 
(q = AM = 0) (Figs. 2a and 3a), as well as the scheme which uses multiple polarized cavity for both the a- dan 
7r-transitions. 

B. Cavity Enhanced Decay 

When the cavity is resonant with the radiative transition from | J, M) to | J — 1, M — q) , the decay rate is enhanced, 
given by the Purcell formula [13] 

r c ,J,M,g = = V^o,.J,M,q (3) 

where T ^j^M,q — TJT^ 2 1S tne ^ ree s P ace decay rate, r\ — the cavity enhancement factor, [ij,M, q the dipole 

transition matrix element, Aj the transition wavelength resonant with the cavity frequency uj c = 2irc/\ c , Q the quality 
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factor and V = J \u mn (r)\ 2 d 3 r the effective cavity volume for the Hermite-Gaussian (m,n) transverse mode function 

We consider a Gaussian mode in a symmetric confocal resonator (aligned along z-axis, see Fig. 1), we have [14] 
the cavity spacing L = (s + \)\ c /2 for longitudinal mode order s, field waist w = \J 4^ = X C ^J (s + \)/^ intensity 

Ko(r)| 2 = ■^j ) e' 2{xl+v2)lw ' 2{z) with w 2 (z) = w 2 (l + (^f) 2 ). From these, the effective volume and the enhancement 
factor is rewritten as 

V = L W 2 n/2 = ±f = ( S + 1 -) 2 \l/16 (4) 
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Large enhancement factor is obtained by using a superconducting cavity with large Q and a low order mode s. The 
free and enhanced decay rates are dependent on the J and M through the numerical factor in the matrix elements. 



C. Shielding thermal photons 

In microwave regime, the mean number of thermal photons n(w c ) is not negligible. It enhances the spontaneous 

2 3 

decay rate as T 0t j t M,q — ^s^nC^t (2n(o; c ) + 1). However, it reduces the cooling efficiency since the decaying states 
cannot be completely depopulated due to simultaneous incoherent excitations of the thermal photons. The usual 
way to reduce the thermal photons is by cooling of the apparatus. The inner surface of the apparatus is coated with 
microwave photon absorbing material like graphite while the outer surface can also be shielded with a photonic band- 
gap structure [15]. Thus, we can set n — and write the free space decay rate for 7r-transition \ J,M) — > | J — 1,M) 
and cr-transition | J, M) — > | J — 1, M =F 1) [16] respectively as 
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r o .j,M, 9 = M)2. 818 7 x 10 46 ( in S.I.) (6) 

where \\fi\\ is the reduced electric dipolc matrix clement for the vibronic state, and g(J,M) = (2j-^)(2j+ij an< ^ 
^t^-i)!^ /+i ) +2 ^ ^ or ^ e ^-transition and cr-transition respectively. 



D. Electric field tuning 

In order to enhance the decay of each transition, we can either tune the cavity or the levels into resonance with 
each other. In principle, we can construct a tunable cavity by changing the distance between the mirrors. First, the 
molecules are put into the cavity tuned to enhance the spontaneous emission of the highest rotational state, from J max 
to J max — 1- Then, the cavity is tuned to the next lower transition and so on until all the populations are brought 
to a single rovibrational state J = 0. Here, L changes from (s + j)Aj max /2 to (s + j)Xq/2. However, this approach is 
not very feasible in practice [17]. 

It is more realistic to tune the transition levels into resonance with cavity frequency using an external electric 
field. For AM — transitions, photons are most likely to be emitted perpendicular to the dipole. To maximize the 
molecular dipole-cavity field coupling strength. and the spontaneous decay rate, the dipole should be parallel to the 
cavity electric field (in x-y plane). Therefore, the tuning electric field is applied along the x-axis so that it aligns the 
dipole along the cavity field (see Fig. 1). Here, the AM = ±1 transitions can be enhanced too as photons are emitted 
in all directions, although the most probable direction is parallel to the dipole. If the tuning electric field is along the 
cavity axis, the AM = transition may not be optimumly enhanced because no photon is emitted exactly along the 
dipolc. However, there is still some enhancement from off-axis modes because the microwave cavity is essentially a 
closed cavity which provides a large solid angle of mode confinement. 

A pair of closely spaced electrostatic plates can be introduced into the cavity to produce a tunable homogeneous 
transverse electric field Ej perpendicular to the cavity axis (say along x-direction) for Stark tuning. Thus, the edges 
of the superconducting cavity mirrors can be shielded from the strong electric field using a metallic casing or coating 
to maintain a constant high Q value. 

The energy separation between states | J, M) and | J — 1, M — q) in the electric field Ej t M,q is 
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AC,i,M,, = hc2J{B e a e (n + -)} + ^™ 9 J (h,M fj-i,M- g ) (7) 

where f J>M - /j_i,m- 9 = ./(j+i^j-i^L+a) _ ( f- 1 ) j*2 j - 3H 2 j+ 1 ) and ^ is tne electric dipole moment. 

From Eq. 7, it is possible to keep a pair of transition at cavity resonance A£" JM(J = by varying the electric 
field for each set of J, M, q. Setting n = 0, the required electric field for each transition can be predicted as 



Ej, M , q = -\l- f ^ {f - J(2B e a e )} (8) 

JJ,M — Jj-l,M-q A c 

We need to fix the cavity parameters through A c such that the argument in the square root is positive. In general, 
f.JM — fj-i,M-q can be positive or negative depending on the transition, so it may not be possible to use a single 
value of A c to tune all transitions to cavity resonance. It is only possible to tunc those transitions with negative values 
of fj t M — fj-i,M-q (the 'negative signs' transitions in Figs. 4a and b) to cavity resonance by choosing the cavity 

parameter A c which gives a negative nominator in Eq. 8, j- < J(2B e — a e ) or L > 2 j(2B -<* ) m or dcr to obtain 



real values of E JM , q - This can be satisfied for J ^ 1 and J > ^ \ (\ + 6M 2 + + 3M 2 + 36M 4 )) if we set 
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(2B e - a e ) or L = , op 2 . (9) 



and Eqs. 6 and 8 become 

r o ,j,M, g = |H| 2 (2B e -a e ) 3 3 (J,M)2.8 1 8 7 x 10 46 ( in S.I.) (10) 



E, M „ = h A l ^r^ ( w) (id 

M V JJ,M - JJ-l,M-q 

Condition Eq. 9 applies to the 7r-transitions (Fig. 4a) and the cr-transitions (Fig. 4b) where /j ; m — fj-i,M-q are 
negative. It also applies to | J = 1, M = 0) — > |0, 0) although / 1)0 — / ,o is positive (Fig. 4a) because Eq. 11 vanishes 
for J = 1. 

On the other hand, the "corner" 7r-transitions \ J,M — ±J) — > \ J — 1,±(J— 1)) have positive denominator fj t ±j 
— /j_i,±(j_i) = (2 7+i)(2j+3)j( 7+1) • ^ nc cay ity value A c defined by Eq. 9 gives negative nominator and cannot be 
used to tune these transitions. The "corner" transitions and those 7r-transitions (Fig. 4a) with 'positive sign' can be 
tuned for cooling by using a different cavity parameter which satisfies a positive nominator in Eq. 8, j- > J(2B e — a e ). 
A reasonable choice is 

^ = J max (2B e - ct e ) (12) 

and we have 

{2B e - a e fg(J, M)2. 8 1 8 7 x 10 46 ( in S.I.) (13) 




{ Jmax - J} (14) 



If the cavity only supports linearly polarized photons, only the parallel or it- transitions (AM = 0) can be enhanced. 
For the moment, we consider only the 7r-transitions, from the states | J, \M\ < O) to | J— 1, \M\ < Q). The transitions 
are tuned to cavity resonance at each time using the electric field, while all other Zeeman states are not in resonance 
with the cavity due to the different Stark shifts. Only those molecules in the states | J, —SI < M <il) are rotationally 
cooled down to a number of 20 + 1 states in the ground level, |fi, — ft < M < ST). For example, with OH molecules 
in state 2 n n= i (Fig. 3a), rotational cooling leads to 2 internal ground states. For spinless molecules like CsF in the 
ground electronic state X 1 Sr 2= o (Fig. 2a) the molecules can be cooled to a single internal state, |0, 0). The number of 
tuning steps required for AM — transitions is (Sl + l)(J max — ST) for integer SI, and (Sl+ |)(J m ax — SI) for half- integer 
SI. The tuning electric fields calculated from Eq. 11 are shown in Figs. 2b and 3b for different J of the decaying 
state. 
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From Figs. 4a and 4b, we see that it is possible establish a scheme to cool the molecules in most states by using a 
multiple field polarizations cavity which can support all the three polarized photons (linear, a~ and er + ) corresponding 
to AM — 0, ±1 transitions respectively. This enables the molecules in all the Zeeman states except the 'corner states' 
| J, ±J), to be Stark-tuned for cooling into a single ground level (or single internal state for case O = 0). Due to the 
quadratic dependence of the Stark shift on M, two pairs of states can be Stark-tuned to cavity resonance at each 
time; but only one pair of states for M = M = transition. 

The total number of steps for cooling with multiple polarizations is generally given by ^(Jmax — Jl)(J m ax + H + 3) 
for integer O and ^(J max — !J)(J max + + 2) for half-integer O while the total number of transitions is J max (^max + 
2) — 0(0 + 2) . The possible sequence of tuning is shown in Fig. 5 for the molecules with spin S = 0. The two 
sequences: A and B correspond to different cavity parameters given by Eq. 9 and Eq. 12 respectively. Sequence A 
requires ^(J max — J!)(J ma x + + 1) steps for integer O and ^(J^ax — ^ 2 ) steps f° r half-integer O while sequence B 
requires ( J max — ^) steps. Sequence A cools J max + + 1 times more states than sequence B but the transitions of 
both sequences are complementary. The tuning electric fields for sequence A and sequence B are calculated from Eq. 
11 and Eq. 14 as shown in Figs. 6a and 6b respectively. 

E. Doppler broadening 

The cavity lincwidth k = uj c /Q must be larger than the maximum Doppler shift in order for the moving molecules 
to stay in resonance with the cavity u> c /Q > lud = w c P max /Mc. For thermal ensemble of OH molecules at about 
10K (70ms -1 ), the maximum quality factor we can use is Q Mc/P max w 10 6 . This value just gives the highest 
enhancement factor in the bad cavity regime. So, it is not helpful to use larger Q. The momentum distribution of 
the molecules is essentially unaffected by the small microwave photon recoil momentum from rotational spontaneous 
emissions unless the molecules are ultracold. 

F. Confinement 

The maximum diameter (radial dimension) for confocal cavity geometry can be estimated as D — \fiL (see Fig. 
1). For a typical microwave cavity dimension of 1mm and molecules with velocity of say 70ms -1 , the interaction time 
of At — w/u — 0.01ms is much shorter than the enhanced lifetime. Therefore, we need to confine the molecules 
within the cavity. 

Figure 1 shows the proposed rotational cooling apparatus with cavity and electric octupoles. Higher order elec- 
trostatic poles can be used to confine the low field seekers with kinetic energy below about IK [18]. This background 
field also prevents the Majorana flop into untrapping state. 

It may be possible to confine molecules using a non-stick solid material like teflon which has low surface adhesion to 
certain molecules. However, this is only applicable to room temperature gas and not translationally cold gas because 
collisions with the solid wall will lead to translational thermalization. Since the molecules in the solid material do not 
rotate and have no angular momentum, collisions with the surface do not cause rotational transition via conservation 
of angular momentum. The main advantage of solid material confinement is of course the infinitely long trapping 
time. 

G. Numerical estimates 

We estimate the cooling time with typical applied tuning field for two types of molecules only for the 7r-transitions 
involving the middle states (—0 < M < O). First, we consider CsF molecules in state (Fig. 2a) with B e = 

0.1844cm" 1 and a e = 1.18 x 10" cm , and the electric dipole moment /i = 7.87Dcbyc [21]. The required electric field 
is in the order of l0 7 Vm~ 1 (Fig. 2b). The cavity spacing from Eq. 9 is L = 2.04cm for mode order s = 1. The cavity 
diameter is D ~ 3.5cm. If the distance between opposite electrostatic poles which produces the homogeneous electric 
field Ej is 1cm, the required voltage is around 100/cV! The free space decay rate is r 0: j, ,o = (2 7-i'|(2j+i) 9- 65x 10 _7 s _1 
and for Q = 10 6 , the enhancement factor of 77 = 0.54 x 10 6 boost the decay rate to r c ,j :0 ,o = 0. 52 ^ 2 j-i)(2J+\) s ^ "• 
Supposing that a duration of tj = 4/r Ci j j o,o ~ 8 ( 2J ~ 1 j( 2J+1 ) j s allocated for each stage of decay from level J, then 

the total duration for enchanced decay from J max = 5 to J Q = is t to t = tj — Yl & ( 2J ~ 1 jf> 2J+1 ) rj 150s (<~ 2.5 

J=J + 1 J=l 
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minutes). Since the trapping volume is macroscopically large, many molecules can be confined and this allows the 
use of sufficiently dilute gas to reduce loss rate of molecules due to state changing collisions. Due to the absence 
of Majorana flop, a significant number of molecules should remain confined for long enough time for the rotational 
cooling to take place. 

Next, we estimate the essential cooling parameters using OH molecules (with the levels in Fig. 3a), which has 
a large rotational constants B e — 18.91cm -1 and a e = 0.724cm~ 1 [19], and the static dipole moment of \i = 
l.6676Debye [20]. From Fig. 3b, we find that the required electric fields are around 10 10 T^m _1 . Even for cavity 
dimension of 0.1mm, the required voltage is very large V ~ lOOOfeF and the free space decay rate is r„ ; j>/,o = 
( (2j-i)(2j+if 4 452 1 x 10~ 2 s _1 - For a small Q = 10 3 and assuming J max = 5, we have r\ = g^np = 540. 93 and 
U , - T AT' 1 - V l 2(2./-i)(2j+i) 6? 

Hot — c,J,M,0 ~ Z^ i ^24(J+l/2)(J-l/2) ~ 

For OH, the required electric field may be too large to be realized in practice and it is due to the unusually large 
rotational constant B e . However, it requires a small cavity Q and the cooling time can be considerably shortened if 
the highest Q in the bad cavity regime is used. 



IV. CONCLUSIONS 



We have proposed a rotational cooling scheme for polar molecules using Stark-tuned internal levels in a supercon- 
ducting microwave cavity. Homogeneous electric field is applied for Stark tuning each transition to cavity resonance 
to enhanced rotational spontaneous emissions. Sequential tuning and the use of the cavity which supports multiple 
field polarizations enables all molecules to be cooled towards the ground rotational level. Numerical estimates for 
CsF molecules show that the scheme requires high voltage which can be realized using the state-of-the-art technology. 
However, molecules with unusually large rotational constants requires electric field beyond the current capability. The 
typical cooling time of one minute with a large Q enhancement in dissipative cavity regime can be realized for most 
molecules with moderate rotational constant and large dipole moment. 
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FIG. 1. Schematic of the apparatus for rotational cooling using cavity enhanced dissipation and electric field for resonance 
tuning and 3-d confinement. 
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b) 

FIG. 2. Spinless (S=0) molecule like CsF in state |X 1 Eo,n = 0} : a)Rotational levels and the Zeeman states with a single 
ground state. Only the parallel transitions 1 to 3 with AM — are shown, b) Electric fields for tuning the CsF molecules from 
the upper state | J max = 10, M = 0) to the lowest state | J = 0, M' — 0) 
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|J=l/2> 



M--1/2 M=I/2 



a) 




b) 

FIG. 3. OH molecule in state \X 1 H 1 /2,n — 0} : a)Rotational levels and the Zeeman states with two ground states. Six parallel 
transitions 1 to 3 are shown, b) Electric fields for tuning the OH molecules from the upper states of |J max = 9.5, M = ±0.5) 
to ground states \J = 0.5, M' = ±0.5) 



J-4) 
|J=3> 



|J=2> 

|J=1) 
J-0> 



M=0 
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|.I-0> 



M=0 

AM-0 transition 
AM= ± 1 transitions 



b) 

FIG. 4. The signs of /j.a/' — fj-i,M in Eq. 8 for: a)7r-transitions (AM = 0) and b)7r-transitions (AM = ±1). Due 
to the quadratic dependence of the Stark shift on M, the transitions | J, \M\ ±1) — » \J — 1, \M\) are symmetrical to 
|J,-|M|-1>-|J-1,-|M|>. 
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|J=3> 




Sequence A : Sequence B : 

ji-transition (AM-0) ► a-transition (AM- +1) 

* a-transition (AM- +1) 

FIG. 5. A scheme of tuning sequence for rotational cooling using a multiple field polarization cavity. Sequence A (transition- 1 
to transition-6) is used with the cavity which satisfies j- = (2B e — a e ). Sequence B is composed of the 'corner transitions' (a 
to c) for cavity with y- = J max (2_B e — a e ). 
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b) 

FIG. 6. Electric fields for tuning the CsF molecules using multiple field polarization cavity. a)Sequence A - The points 
corresponding to each J are displaced slightly from from each other for clarity and connected by line to show the tuning 
sequence which starts from the largest electric field value. The points 'o' correspond to the ir transitions. b)Sequence B - The 
'o' line is for J max = 4 while the '*' line is for J max = 10. The abscissa J labels the upper rotational state for each transition. 
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